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Abstract 

The traditional C-V model had a problem of intensity 
inhomogeneous in segmentation for images. This paper 
proposed a modified C-V model based on local gray 
information, called Local Gray Chan-Vese (LGCV) model, 
through using local gray image information, this model 
could segment inhomogeneous image in less number of 
iterations. Experimental results also prove the effectiveness 
of this method. 
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Introduction 

Ultrasound image segmentation was a typical problem 
in medical image segmentation due to the poor 
contrast and gray uneven which were caused by the 
speckle noise. After decades of development, a lot of 
methods have been proposed for image segmentation 
technology, such as, segmentation method based on 
Region, segmentation method based on Edge, Region 
segmentation method combined with edge, 
segmentation method based on Graph cut theory and 
so on. 

At present, geometric active contour model has been 
widely application. Based on the MS model, Chan- 
Vese proposed C-V model, which was a segmentation 
method based on the regional statistical information 
and without depended on the gradient information of 
the image and with the advantages of not sensitive to 
light and noise. 

According to the segmentation problems of the gray 
level uneven image in C-V model,this paper increased 
local gray level information energy item based on the 
C-V global energy model.lt could better adapt to the 
characteristics of the uniform and nonuniform gray 
level image by adjusting the coefficient of the global 
energy item and the local energy item. 


The Improved C-V Model 

The Local Item and the Regularization Item 

For the segmentation of gray level nonuniform image. 
This paper presented the integration into the local 
gray level statistical information in the C-V model, 
which contains the LGCV model of the local 
energy E L . 

In the images with intensity inhomogeneity, the 
uneven grayscale effect was slowly changing, so the 
gray level kept the uniformity characteristics in a small 
local region.In this paper, it chose the mean of local 
gray level distribution as a statistical information, then 
made difference between the image convolution and 
the original image, which enhanced the intensity 
contrast between the target and the background 
greatly. The regional grayscale average value could be 
replaced by the local grayscale average value through 
combining with level set method, and obtained the 
mathematical expression of E L , as follows: 

E L id x ,d 2 ,C) = f inside(C) | g k *u 0 (x, y) - u 0 (x, y) - d x | 2 dxdy 

+ L, side(C) \Sk * «o (*. y) - «o (*. y) - d 2 1 2 dxdy 

( 1 - 1 ) 

Where g* was the average convolution operator with 
size window of k , d x and d 2 were the gray level 
average value in internal and external curve of the 
difference (g* * u 0 (x, y) - u 0 (x, y)) between the image 
convolution and the original image,respectively. 

When the evolution curve approximated the real 
target edge of the difference image in time, the value 
of the local energy item would decline rapidly and the 
solution after the local energy reached minimum was 
the target contour C* . 

inf {F X L (C) + F 2 l (O) - 0 kF x l (C * ) + F 2 L (C* ) (i_ 2 ) 

Where F { L (C) represented the first item in equation 
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(1-1), F 2 (C) represented the second item in equation 
( 1 - 1 ). 

It introduced the level set function $(x,y) to equation 
(1-1), so the local item e l could be rewrote as follows: 

E l ( d } ,d 2 ,<f>) = \ n | (g k * Uq (x, y) - U q (. x , y)) -d x | 2 H (</>(x, y))dxdy 

+ In \g k * «o (*, y ) “ «o (*> >0 - d 2 1 2 C 1 - H y)))dxdy 

(1-3) 

Furthermore, the length penalty item L(C) was added 

to E l regularization item expression for keeping the 
smooth of zero level set and avoiding to produce the 
small and isolated area in the final segmentation 
result, which also made the curve keeping as short as 
possible in the evolutionary process.Assumed C was 
a smooth and closed plane curve: C(p ) :[0,1] -► Cl , 
parameter p £ [0,1], the length of curve C could be 
expressed as follows: 

L(C) = \$\ c dp (1-4) 

Using the zero level set of level set function instead of 
the curve C , thus L(C) was expressed as follows: 

L((/> = 0) = In | V// (</>(x, y))\dxdy = \ n 8(<p(x,y))\V <p(x,y)\lxdy 

(1-5) 

In addition, an energy penalty item was added to the 
regularization item expression for keeping the level set 
function approximated the signed distance function. 
The item was expressed as follows: 

m = (1-6) 


fidelity item in C-V model, as follows: 

E g (c, , C 2 , 4>) = J n \u 0 ( x , y) - c, | 2 H(</>(x, y))dxdy + 

2 ( 1_ °) 
l n \u 0 (x,y)-c 2 \ H(<t>(x,y)))dxdy 

The second part was the local gray level item showed 
in (1-3), it was necessary to choose a appropriate size 
window k for the average convolution operator g k . 

The fidelity item could induce the evolving curve to be 
close to the target edge quickly, and the local item 
could constraint the curve to stop at the real edge 
contour of the object, so the local item should be used 
combining with the global item. The last, added the 
regularization item to the model and the total energy 
function of LGCV model was expressed as follows: 

E (c x , c 2 , d l , d 2 , (fk) — cc • E (c x , c 2 , (j>) ^ ^ 

+/3-E L (d l ,d 2 ,</)) + E R (</)) 

Where a and (3 were two positive parameters for 
measuring the effect of the global item and the local 
item. In the actual application process, the value of a 
and P should be set according to the gray level 
distribution of image: (1) For the image with uniform 
gray level distribution, the value of a should 
approximate to p , which meaned that the global item 
and the local item had the same effect; (2) For the 
image with nonuniform gray level distribution, the 
value of a should greater than p , which meaned that 
the effect of the local item was strengthened. 

Expansion the equation (1-9) to equation (1-10), as 
follows: 


In fact, the energy penalty item P{(P} amounted to a 
measurement to measure the similarity between level 
set function ^ and signed distance function, which 
played a key role in avoiding reinitialized process.And 
then, obtained the regularization item E R which was 
expressed as follows: 

E R {p) = p- L{(j) = 0) + Pip) 

= M ■ Jn S{(j>(x, y)) | V <j>(x, y)\ dxdy + J fi 1(| ' V <j>(x, j>)| - 1) 2 dxdy 

(1-7) 

Where // was a parameter of controlling the length 
penalty item,if // was small relatively,the LGCV model 
could segment the smaller area object,and vice versa. 

LGCV Model Segmentation Algorithm 

LGCV model consisted three parts mainly, the first 
part was the global energy item which composed of 


E IGCr (qxiri x ,d 2 ,t) 

= fn (« • h (*, y) - 5 f + P ■ |g* * «0 (*> y) -% W l) - d \ f )H(i/>(x, y))dxdy 

+ J n (« ■ K (*> y) c i | 2 + P • |& * M o (y y) - % y) - d i ( X 1 - H (^( x > y))) dxd y 

+ (,M ■ f Q 8{<Hx, y)) \V<j){x, y) \ dxdy + f Q j (|' V <f>(x, y) \ - 1) 2 dxdy ) 

( 1 - 10 ) 


Where c x (<j>) t , d x (</>) and d 2 (<f>) were expressed 
respectively as follows: 


= j n « 0 (x,y)//(^(x,y))dxdy 
' In H(</>(x,y))dxdy 

= In «o (v l)(i - riQ/Ei 

In(l -H(</>(x,y)))dxdy 

d ( = L (gk * “o ( x > y) z M o ( x ’ yypEl 

1 \ (1 H(<t>(x,y))dxdy 

d = In (gk * «o y) ~ M o (V jOXl ~ H(tf>(x, y)J)dxdy_ 

In (1 - ri(<f>(x, y)))dxdy 


( 1 - 11 ) 

( 1 - 12 ) 

(1-13) 

(1-14) 
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In the light of the above value obtained of c lf c 2 , d x 
and d 2 ,the level set evolving equation of energy 
functional e lgcv could be expressed as follows: 


Q ^ (h 2 9 

— = -c,) + /3-(g k *u 0 (x,y)-u 0 (x,y)-d l ) ) 

ot 

+ («-(«o - c if + P ' (Sk *u 0 (x,y)-u 0 (x,y)-d 2 ) 2 )] 

+ [ju ■ 8(<t>)div(^-) + (V V - div { ^-M 


The Experimental Results and Analysis 

In the LGCV model, (3 - 1 ,and dynamically adjusted 
the value of the global item parameter a according to 
the actual situation of the image. In the experiment, 
when a < 1 ,it adjusted all the gray level nonuniform 
images,when a= 1 ,it adjusted all the gray level 
uniform images.As similar to a ,the length penalty 
item parameter was adjustable.In this experiment, 
ju= 0.01*255 2 corresponded to the small object 
segmentation, and ju= 0.1*255 2 corresponded to the 
large object segmentation,so it only adjusted a and ju 
in the following tests. 

Figure 2-1 showed the segmentation results of C-V 
model and LGCV model to an artificial simulation 
image respectively. Parameters value in Figure 2-1: 
a = 0A , J3 = 1 (LGCV model), ^=^=1 , 

ju= 0.01x255 2 (C-V model and LGCV model). 


(a)original image (b)initial contour 


(c)C-V algorithm (d)LGCV algorithm 

FIG 2-1 SEGMENTATION RESULTS OF SYNTHETIC IMAGE 

Figure 2-2 were the medical ultrasound images. For 
saving operation time, the region of interest in the 
images should be cut out as the figure 2-3(a) and 2-4(a), 
and then performed segmentation operation. 

Added the energy penalty item to the C-V model in 
order to the fair comparison, so the C-V model had no 






need to reinitialize in the evolutionary process. 
Parameters value in Figure 2-3: a = 0.5, /3 = \ (LGCV 

model), 4=^=1, /u = 0.01 x255 2 (C-V model and 
LGCV model). 



(a) liver cancer (b) breast tumor 


FIG 2-2 CLINICAL ULTRASOUND IMAGES 



(a) original image (b) initial image 



(c) C-V algorithm (d) LGCV algorithm 


FIG 2-3 SEGMENTATION RESULTS OF LIVER CANCER IMAGE 
BASED ON CV MODEL AND LGCV MODEL 

Parameters value in Figure 2-4: a = 0.8, J3 = \ (LGCV 

model), ^=^=1 , /j = 0.1x255 2 (C-V model and 

LGCV model). 



(a)original image (b)initial image 



(c)C-V algorithm (d)LGCV algorithm 

FIG 2-4 SEGMENTATION RESULTS OF BREAST TUMOR IMAGE 
BASED ON C-V MODEL AND LGCV MODEL 


Conclusions 

Based on the C-V model, according to the image 
segmentation method of using the global information 
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in C-V model, and regarded the fidelity item of C-V 
model as the global energy item in the new model, 
then increased the local gray level information, the 
local energy item and the regularization item, and 
finally proposed the improved model, LGCV model. 
The model could obtain a good segmentation result in 
gray level uneven image by adjusting the weight of 
the global item, the local item and the regularization 
item. Compared to C-V model, the experimental 
results show that LGCV model has a better 
segmentation effect in the emulation image and the 
actual medical ultrasound image with uneven gray 
level. 

ACKNOWLEDGMENT 

The work was supported by the Foundation of 
Guangdong Sci-Tech Plan. 

REFERENCES 

Chan TF,Vese LA. " Active contours without edges." IEEE 
Transactions on Image Processing, 10(2001):266-267. 
D.Ziou, S.Tabbone. "Edge Detection Techniques-An 
Overview." Department of Math Information. University 
of Sherbrook, 1997. 

Kass M, Witkin A, Terzopoulos D."Snakes: active contour 
models." International Journal of Computer Vision, 
1(1987): 321-331. 

Mumford D,Shah J. "Optimal approximation by piecewise 
smooths functions and associated variational problems." 
Communications on Pure and Applied Mathematics, 
42(1989):677-685. 


Osher S, Sethian J A. F"ronts propagating with curvature- 
depen-dent speed: Algorithms based on Hamilton-Jacobi 
formulations." Journal of Computational Physics, 
(79)1998:12-49. 

P.K.Sehoo, S.Soltani, A.K.C.Wang, etc. A Survey of 
Thresholding Techniques. Computer Vision, Graphics 
and Image Processing. 1988.41(2):233-260. 

Minfen Shen Shantou, Guangdong, 1957.9. He received the 
M. Eng.degree in signal processing and Ph. D. degree in 
communication and information systems from the 
University of Science and technology of China, Hefei, China. 

After graduation, he joined the Department of Electrical and 
electronic Engineering, Shantou University, Shantou, China, 
where he is currently a Professor.He is also an adjunct 
professor with the Department of Science and Technology, 
USTC.His research interests include random signal 
processing,computational intelligence, biomedical signal 
analysis,nonlinear signal and image processing, and chaos. 

Lixia Zhang Xiaogan,Hubei,1991.11. She is a graduate 
student in Shantou University,China. Her research interests 
include medical signal and image analysis. 

Lisha Sun Wuhan, Hubei, 1958.7. She received the B. 
Eng.degree in Electronic Engineering from the University of 
Geology Science of China, Wuhan, China. 

After graduation, she joined the Department of Electrical 
and electronic Engineering, Shantou University, Shantou, 
China, where she is currently a Professor. Her research 
interests include medical signal and image analysis. 

Zhifei Su Jingzhou ,Hubei, 1987.8. He is a graduate student 
in Shantou University, China. His research interests include 
medical signal and image analysis. 


34 



